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Measurement of parity-violating spin asymmetries in production at midrapidity 

in longitudinally polarized p+p collisions 
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We present midrapidity measurements from the PHENIX experiment of large parity-violating 
single spin asymmetries of high transverse momentum electrons and positrons from W ± /Z de¬ 
cays, produced in longitudinally polarized p+p collisions at center of mass energies of -^8=500 and 
510 GeV. These asymmetries allow direct access to the anti-quark polarized parton distribution 
functions due to the parity-violating nature of the VE-boson coupling to quarks and anti-quarks. 

The results presented are based on data collected in 2011, 2012, and 2013 with an integrated lu¬ 
minosity of 240 pb _1 , which exceeds previous PHENIX published results by a factor of more than 
27. These high Q 2 data probe the parton structure of the proton at W mass scale and provide an 
important addition to our understanding of the anti-quark parton helicity distribution functions at 
an intermediate Bjorken x value of roughly Mw/y/s = 0.16. 

PACS numbers: 14.20.Dh, 25.40.Ep, 13.85.Qk, 13.88.+e 


The determination of the contributions of partons to 
the spin of the proton has inspired significant theoret¬ 
ical and experimental effort for several decades [M3|. 
The quark contribution to the nucleon spin has been 
deduced through measurements in polarized inclusive 
deep-inelastic scattering (DIS) and semi-inclusive deep- 
inelastic scattering (SIDIS) experiments pi fT^UT5| . Al¬ 
though the overall quark contribution (AS = A q + A q) 
has been well-determined through DIS experiments (in 
the range 10~ 3 < x < 1), the contributions from sea 
quarks separated by flavor (determined through SIDIS 
experiments) are comparatively poorly known. Data 
from HERMES and COMPASS pi [Tf3] provide con¬ 
straints on the contribution from the sea quarks, however, 
uncertainties in fragmentation functions and the low en¬ 
ergy scales of fixed target experiments limit the accuracy 
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with which these measurements can quantitatively deter¬ 
mine the sea quark contribution m- As such, an inde¬ 
pendent measurement using a different technique (1(5] to 
determine the contribution from different flavors of sea 
quarks is desirable. 

The use of VE-boson production provides just such a so¬ 
lution. Parity is maximally-violated in the W couplings 
to quarks and leptons, so W^ 1 production in p+p col¬ 
lisions proceeds only by coupling to left-handed quarks 
and right-handed anti-quarks ( ulAr —I W + and (IlUr —I 
W~). By measuring decay leptons in the final state, the 
flavor and helicity state of the colliding quarks can be 
determined [MU- Asymmetries measured in W ± by 
reversing the helicity of a colliding proton are sensitive 
to the individual quark/anti-quark helicity parton distri¬ 
bution functions (PDFs) (A u, Ad, A u and Ad). More¬ 
over, the energy scale for these events, of the order of the 
VE-boson mass, allows for small and precisely calculable 
higher-order corrections. 

We present results for the parity-violating single-spin 
asymmetry Ar for p + p —i W ± /Z + X —> + X' 
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at midrapidity from 2011 2013 PHENIX data at the 
Relativistic Heavy Ion Collider (RHIC). These results 
relate to an intermediate Bjorken x value of roughly 
Mw/yfs = 0.16. Initial measurements at RHIC in 2009 
accumulated 8.6 pb” 1 by PHENIX |JJj and 12 pb” 1 by 
STAR [lOl HI] . Here, the total integrated luminosity is 
240 pb" 1 at = 500 GeV in 2011, and at 510 GeV in 
2012 and 2013 [22]. Proton-beam polarizations were also 
considerably improved from ^0.39 in 2009 to 0.50-0.56 
in 2011-2013. 

The measurements are performed with the two 
PHENIX central arm spectrometers. Each arm covers 
|A0| = 7t/2 in azimuth and |ry| < 0.35 in pseudorapid¬ 
ity. A comprehensive description of the PHENIX detec¬ 
tor at RHIC can be found in [23]. The major detector 
subsystems used for this analysis are the electromagnetic 
calorimeter (EMCal) and the drift chamber/pad cham¬ 
ber tracking system. Two beam-beam counters located 
at ± 144 cm from the collision point along the beam line 
and covering 3.1 < |? 7 | < 3.9 were used to define the min¬ 
imum bias trigger and to measure the relative luminosity 
between different colliding bunch pairs. 

The data were collected with an EMCal-based trig¬ 
ger [23] with nominal energy threshold of 5.6 GeV, 
which was fully efficient for e± with transverse momen¬ 
tum pj, >10 GeV/c. The p? was determined from the 
energy deposited in the E MCal with energy resolution 
t te/E = 8.1 %/^/E{GeV) © 4%. The energy resolution 
was determined from the pr-dependence of the widths of 
reconstructed ^ r° and rj meson mass peaks. The same 
7 r° and p meson mass peaks were used in the energy 
calibration of the EMCal, and were continuously moni¬ 
tored. Similar to our previous analysis [5] and test beam 
data results [23;, the EMCal energy scale was confirmed 
to within 2.5%, for the energy range analyzed with this 
data. A loose time-of-fhght cut was applied in the anal¬ 
ysis to remove noncollision background. 

The tracking system was used for collision vertex re¬ 
construction, track charge sign determination, and back¬ 
ground suppression. The main tracking detector, the 
drift chamber (DC), spanning the radial distance 2.02- 
2.46 m from the beam line, measured the charged track 
bending in the axial magnetic field of the PHENIX cen¬ 
tral magnet, with a field integral of 1.15 Tm. The 21 - 
coordinate for the tracks was obtained from the pad 
chambers situated behind the DC. Reconstructed tracks 
were matched with high energy clusters in the EMCal 
within a cone angle of 0.02, retaining >99% of real e ± 
tracks, as determined from simulations. The coordinate 
information from both the calorimeter and the tracking 
system was used to determine the 2 :-vertex of the event, 
and only events with | 2 r| <30 cm were used in the anal¬ 
ysis. 

The charge sign of a track was determined from 
the bending angle odc , which is inversely propor¬ 
tional to the track transverse momentum (aDc( mr a'd) = 
92 /p T GeV/c). A region corresponding to |aDcl < 1 
mrad was removed in order to minimize the possibility 


of charge misidentification. This led to <3% loss of e ± 
from IV-boson decays. To further eliminate the charge 
sign ambiguity in the DC track reconstruction, the re¬ 
gions in the vicinity of anode wires were removed from 
the analysis, reducing the DC acceptance by ~15%. The 
remaining opposite charge contribution to the W~ (W + ) 
signal was 2% (0.4%), as determined using the DC reso¬ 
lution of 1.4 mrad and bdc convoluted over the W decay 
e ± pt distribution. The result is consistent with a full 
detector simulation. 

Accurate momentum reconstruction in the tracking 
system requires the precise determination of the beam 
position in the plane orthogonal to the beam line. This 
was measured and monitored using straight tracks from 
special runs with the magnetic field off throughout the 
data taking period. 



FIG. 1. (color online). (Upper panels) Spectra for (a) e + and 
(b) e~ using the EMCal for momentum determination from 
p+p collisions at yA = 510 GeV from 2013. From top to bot¬ 
tom, the curves are: solid [red] is the sum between background 
and signal; shaded band [orange] is the background estima¬ 
tion with the uncertainty from the GPR method; dashed [dark 
blue] is the W ± /Z —>■ e ± signal obtained from simulation nor¬ 
malized to the data; and solid [light blue] is the contribution 
from Z —> e + e”. (Lower panels) Point-by-point comparison 
of the data and the fit result: (datai-fit)/ai, where Ui is sta¬ 
tistical uncertainty of the i-th data point. 

An isolation cut was very efficient at suppressing back¬ 
ground events with a high degree of activity around a 
candidate electron (as would happen for jet events). The 
cut parameter r- lso was defined as r- lso = (££)) / E e , where 
Ei is the j-tli EMCal cluster energy and track pt around 
the electron candidate in a cone with a radius in p and 
t f> of 0.4, and E e is the energy of electron candidate. A 
candidate was kept for the analysis if r; so <0.1. 

Figure [I] shows the resulting yield of electron and 
positron candidates for the 2013 data set. A Jacobian 
peak around pj, =40 GeV/c corresponds to e ± from W 
and Z boson decays. The isolation cut removed about 
90% of the background (as was evaluated from the back¬ 
ground dominated region between 10 and 20 GeV/c) and 
left more than 90% of the signal in the Jacobian peak re- 
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gion untouched (as evaluated from simulations explained 
below). Similar results were obtained for the 2011 and 
2012 data sets. Above 30 GeV/c the remaining candidate 
events after the isolation cut are dominated by W and Z 
decay to electrons/positrons, and by background events 
below 25 GeV/c. This background consists mainly of 
high momentum electrons/positrons from conversion of 
7r°/?7 decay photons, charged pions/kaons, b,c —> e de¬ 
cays and accidental matching between high energy EM- 
Cal clusters and tracks in the DC. The Z-boson contribu¬ 
tion in the signal region above 30 GeV/c was estimated to 
be 7% (25%) for the positrons (electrons) after all anal¬ 
ysis cuts were applied, as determined from simulations. 
The asymmetry of the Z has been estimated theoretically 
using the DSSV08 PDF sets and measured by the STAR 
collaboration [IT] to be —0.07±0.14. 

Experimentally, the longitudinal single-spin asymme¬ 
try is defined as: 


L P N+ + RN~ ’ v ' 

where P is the beam polarization, N + (N ~) is the num¬ 
ber of events in the signal region for the positive (neg¬ 
ative) beam helicity and R is the luminosity ratio (rel¬ 
ative luminosity) between positive and negative helicity 
bunches measured using the minimum bias trigger de¬ 
fined by a coincidence of the two beam-beam counters. 
The relative luminosity between different helicity combi¬ 
nations did not differ from unity by more than 2%. The 
asymmetry calculation was performed for events in the 
Pt range from 30 to 50 GeV/c, which defined the sig¬ 
nal region in this analysis. This range was selected to 
optimize the signal to background. Less than 1% of the 
signal is expected above 50 GeV/c. Asymmetries ob¬ 
tained in this fashion must be corrected for background 
events, which are parity-conserving, in the signal region. 
This dilution factor can be defined as ( A — B)/A , where 
A ( B) is the number of all (background) events in the 
signal region 30 < px < 50 GeV/c. The final asymme¬ 
try values can be obtained by dividing the result by the 
dilution factor. 

The background in the signal region was estimated 
using the Gaussian process regression (GPR) tech¬ 
nique [251(281 to extrapolate the background shape 
from the background-dominated region to the signal- 
dominated region. The major advantage of this method 
is that it does not require an a priori known functional 
form to test against data. At its core, this method al¬ 
lows for the determination of the shape of a set of data 
points with statistical uncertainties using only the data 
themselves. Furthermore, the predictions made using 
this method have a mathematically well-defined Gaus¬ 
sian uncertainty. 

Through our use of the radial basis function (RBF) 
kernel [25] [26], we assume a smooth (infinitely differen¬ 
tiable) shape for the background. The background shape 
was constrained from data points in the px ranges 10-22 
GeV/c and 60-65 GeV/c, where the signal contribution 


is expected to be negligible. Although bins in the range 
60-65 GeV/c don’t contain any events, they still improve 
the precision of the background evaluation. These empty 
bins were assigned a statistical uncertainty of 1 count. 
The background in the signal region is assumed to vary 
on px scales equal or larger than those in the background 
dominated regions. 

The RBF kernel contains a characteristic length pa¬ 
rameter that is an indicator of how far away from data the 
background extrapolations can be made. For obtained 
characteristic lengths larger than 30 GeV/c, we con¬ 
cluded that our background estimation (based on data 
between 10 to 22 GeV/c and 60 to 65 GeV/c) in the sig¬ 
nal region (30 to 50 GeV/c) has an appropriate statistical 
uncertainty. 

Table [I] summarizes the background contributions with 
statistical uncertainties obtained using the GPR ap¬ 
proach along with the counts in the signal region for each 
data set. The GPR analysis was performed for differ¬ 
ent px ranges for the background estimation and includ¬ 
ing/excluding the constraint between 60 and 65 GeV/c. 
The results were within the statistical uncertainty of the 
full GPR analysis so no additional systematic was added. 

In Fig. [I] the background and signal shapes were used 
to describe the data points. The only fit parameter was 
the normalization for the signal shape. The signal shape 
was obtained from a PYTHIA simulation [29] of W ± and 
Z-boson decays to electrons/positrons, followed by a full 
GEANT3-based m detector response simulation. The 
simulated events were analyzed using the analysis pack¬ 
age used for the data. The fit quality of the data-driven 
background shape plus the simulated signal shape is rea¬ 
sonable for both e~ and e + spectra. 

As a cross check of the background determination, a fit 
to the data using a phenomenologically-motivated modi¬ 
fied power law function as the background shape was also 
performed (^f (pt) = 1 /p^ +l3lnpT ^. The discrepancy be¬ 
tween the central value results from two methods was 
assigned as a systematic uncertainty for the background 
determination. Another source of uncertainty may come 
from the possible systematic discrepancy between the 
data points and the fit result in some px regions (e.g. 
data excess over the fit in the vicinity of px = 30 GeV/c 
in Fig. [Tj) . Following a conservative approach for un¬ 
certainty evaluation, the sum of the signed differences 
between data points and the fit results within the signal 
region was assigned as an additional systematic uncer¬ 
tainty. The final systematic uncertainty was obtained by 
adding in quadrature the systematic uncertainty from the 
two sources discussed above. Using the GPR-estimated 
background contamination in the signal region, the dilu¬ 
tion factor for each data set was calculated and is pre¬ 
sented in Table [H 

The asymmetry calculation was done following two in¬ 
dependent methods. First the asymmetry was calculated 
separately for each polarized beam using Eq. [l] with the 
polarization for the other beam averaged to zero. The 
final result is a weighted average of asymmetries from 
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TABLE I. Number of events recorded for e + and e for 30 < pr < 50 GeV/c and the background contributions, dilution 
factors, and two-beam polarizations for each analyzed data set. 


Lepton 

Year 

Counts 

Background 

Dilution 

Polarization 




Contribution 

Factor 

B 

Y 

e + 

2011 

70 

2.3 ± 2.3 (stat) ± 0.6 (syst) 

0.97 ± 0.03 (stat) ± 0.01 (syst) 

0.51 ± 0.02 

0.50 ± 0.02 


2012 

105 

2.5 ± 2.5 (stat) j ^'5 (syst) 

0.98 ± 0.02 (stat) (syst) 

0.55 ± 0.02 

0.57 ± 0.02 


2013 

669 

18.6 ± 7.3 (stat) ± 14.9 (syst) 

0.97 ± 0.01 (stat) ± 0.02 (syst) 

0.55 ± 0.02 

0.55 ± 0.02 

e~ 

2011 

27 

1.7 ± 1.6 (stat) ± 0.7 (syst) 

0.94 ± 0.06 (stat) ± 0.02 (syst) 

0.51 ± 0.02 

0.50 ± 0.02 


2012 

47 

5.5 ± 4.7 (stat) ± 2.2 (syst) 

0.88 ± 0.10 (stat) ± 0.05 (syst) 

0.55 ± 0.02 

0.57 ± 0.02 


2013 

233 

13.9 ± 5.6 (stat) (syst) 

0.94 ± 0.02 (stat) 1^9 (syst) 

0.55 ± 0.02 

0.55 ± 0.02 


two beams. A likelihood method was also used in order 
to deal with the lower statistics, particularly in the 2011 
and 2012 data sets. 

The two rings at RHIC with counter-propagating 
beams are designated yellow (y,Y) and blue (b,B). The 
number of expected counts p, y t> for the data sample can 
be expressed as: 

H y b = RybN (l+b-A L P B +y-A L P Y +b-y-A LL P B P Y ) ( 2 ) 

where R y b is the relative luminosity between the colliding 
beam helicity configurations, y ( b ) denotes the helicity of 
the two colliding beams and takes the value of +1 (—1) 
for positive (negative) helicity, the parameter N is an 
average count, Pg and Py are the polarizations of the 
two beams, Agg is the double spin asymmetry. The spin 
asymmetries were calculated by maximizing a likelihood 
function defined using Poisson statistics as: 

£= II R (VybjNyb) j (3) 

y=±l,b=±l 

where N y t, is the spin sorted yield. To calculate the 2013 
positive and negative rj bin asymmetries a generalized 
form for these equations was used. 

Table |TT| summarizes the Ag results. Both of the asym¬ 
metry calculation methods employed gave consistent re¬ 
sults for all the data sets. The systematic uncertainties 
were obtained by propagating the systematic uncertain¬ 
ties of the dilution factors to the final asymmetry values. 
A scale uncertainty of 3.5% from the RHIC beam po¬ 
larization measurements is not included in Table HI The 
asymmetry in the background region was also measured 
and for all cases the asymmetry was consistent with zero, 
within uncertainties. 

These results are shown in Fig. [2] with two theoretical 
calculations: (CHE) [21] for the NNPDFpoll.l [FT] and a 
recent calculation [3Tj using the DSSV 14 PDF sets [52] , 
While the DSSV 14 curve was obtained from a global 
fit of DIS and SIDIS data (including recent COMPASS 
results [MM]), the NNPDFpoll.l uncertainty band con¬ 
tains the 2012 STAR El result for flavor separation in 



FIG. 2. (color online). Asymmetry results from the com¬ 
bined 2011 and 2012 data sets for \rj\ < 0.35 (black circles) 
and the 2013 data (red squares) separated into two equal r/ 
bins between -0.35 and 0.35. The green line and shaded re¬ 
gion shows a theoretical calculation using CHE m with the 
NNPDFpoll.l PDF sets [2], while the dashed magenta line 
shows the DSSV14 calculation 1311 . 


addition to DIS data. The theoretical asymmetry calcu¬ 
lations agree with the data within 1.5 a uncertainty of 
the data points. These results will be used to further 
constrain the quark and anti-quark polarized parton dis¬ 
tributions functions at an intermediate Bjorken x value 
of roughly M w /yfs = 0.16. 

Figure [3] shows the combined asymmetry for all of the 
PHENIX data sets and published data from STAR jTTj . 
The two data sets cannot be compared directly, because 
PHENIX measures the asymmetry from W ± + Z decays, 
while the STAR result is solely from W ± decays. The 
comparison can be made through the curves, which ac¬ 
count for the specifics of each measurement. Qualita¬ 
tively, both data sets show the same trend with data 
points below (above) the central value of the theoretical 
prediction for W + (W~ ), for |ry| < 0.5. The W~ dif- 
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TABLE II. Longitudinal single-spin asymmetries, Al, for 
the 2011 and 2012 data sets (combined) spanning the entire rj 
range of PHENIX (\r]\ < 0.35), for the 2013 data set separated 
into two 77 bins, and for the combined 2011-2013 data sets. 

Lepton Data Set (r/) Al 

e+ 2011+2012 0 -0.27 ± 0.10 (stat) ± 0.01 (syst) 

2013 ?? > 0 0.17 -0.38 ± 0.07 (stat) ±0.01 (syst) 

2013 77 < 0 -0.17 -0.35 ± 0.07 (stat) ±0.01 (syst) 

2011-2013 all 0 -0.35 ± 0.04 (stat) ±0.01 (syst) 

e" 2011+2012 0 0.28 ± 0.16 (stat) ± 0.02 (syst) 

2013 ?? > 0 0.17 0.10 ± 0.13 (stat) tg'.oi (syst) 

2013 7 ? < 0 -0.17 0.17 ± 0.12 (stat) (syst) 

2011-2013 all 0 0.17 ± 0.08 (stat) ±0.02 (syst) 



-1 -0.5 0 0.5 1 -1 -0.5 0 0.5 1 

% 


FIG. 3. (color online). Asymmetry results from the com¬ 
bined 2011 -2013 data sets from PHENIX [red] circles and the 
STAR 2011-2012 El W results [blue] stars and their respec¬ 
tive DSSV 14 theoretical predictions. 

ference corresponds to a larger A u contribution in the 
covered x ~ 0.16 range, when compared with the central 
value of the DSSV 14 PDF fit calculation. 

In summary, for high pr e~ and e + from W and Z 
boson decays, PHENIX measured the single spin asym¬ 
metries with more than 27 times higher statistics and 
better polarization compared to 2009 [SJ. These new 


results and the STAR data m will help constrain the 
anti-quark PDFs in a global analysis. Asymmetries cal¬ 
culated from global fits based on previous measurements, 
such as DSSV14 and NNPDFpoll.l, are consistent with 
our data. The use of the electroweak interaction provides 
an independent tool to extract quark and anti-quark he- 
licity contribution. The data presented here are com¬ 
plementary to previous SIDIS measurements and bring 
the field one step closer to elucidation of the proton-spin 
puzzle |T]. 
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